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ABSTRACT 



This thesis exeirdnes the theoretical performance of a 
solid propellant rocket motor which was developed for launching 
small amateur research rockets. The theoretical results are pre- 
sented in the form of two limiting cases concerning the behaviour 
of the two-phase exhaust flow. Actual testins of the motor is 
performed utilizing a specially designed test ris in order to 
compare the results. As wellr optimization of the motor's per- 
formance is investigated . 

The theoretical performance is found to be in stood agree- 
ment with the test results? providing a basis for future desiSn 
of larder endines. 

No sidnificant improvement in overall motor performance 
as a result of modifications is foreseen. 
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CHAPTER 1 



INTRODUCTION 



This thesis has three primary objectives* to theoretic- 
ally analyze the operation of a small solid propellant rocket 
motorf secondlyf to conduct actual testins with which to compare 
the theoretical result? and thirdly» to modify the propellant and 
motor desisfn for optimization of motor performance* 

Secondary objectives are that knowledge of the actual 
impulse (thrust-time traits) of the motor can be instrumental 
in the determination of such parameters as acceleration* velocity 
and altitude of the rocket durins actual fliahts* Finally* the 
research results aenerated by this study would be of value to 
others interested in the field of amateur rocketry* 

1*1 ROCKET MOTOR 

The rocket motor under examination in this research 
study is one of a number of rocket motors developed < largely 
throudh experimentation) durina twelve years of personal amateur 
rocketry * This motor has been used ouite extensively for racket 
fliahts* capable of propellina a 2 ka** 9 cm* diameter rocket to 
an altitude of 500 metres* with a hiah dearee of reliability <95 
percent success rate on a data base of over 50 firinas)* 

1*2 PROPELLANT 

The propellant utilized in the rocket motor is one 
that is employed exclusively by amateur rocket enthusiasts* It 
is not a hiah performance propellant and is not used in prof- 
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essionallu desisned rockets* Therefore» there wss no known data 
available on the performance characteristics of this propellant* 
This necessitated a theoretical 3nal>:<sis coupled with experiment- 
ation on combustion product anal\:isisr combustion temperature 
measurements* burn rate measurement and the effects of varyind 
the oxidizer-fuel (0/F) ratio* 

1.3 EXPERIMENTATION 

Testind of the actual motor was conducted on a specially 
desidned static testind stand* The static testind stand per- 
mitted the motor thrust to be recorded continuously throudhout 
the firind* accomplished by convertind the thrust to an electric 
analodue sidnal which was amplified* diditized by an A/D converter* 
then stored in a computer for processind* 

Optimization of the propellant involved testind the 
effects upon its performance by varyind 0/F ratios* The motor 
optimization consisted of modification to the nozzle desidn* in 
order to attempt to increase the nozzle efficiency* 





Fiisure l.l* Rocket launch eauipped with small solid motor 
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CHAPTER 2t THEORY 

2«1 Solid Propellant Rocket 

A typical solid propellant motor (fidure 2»1) is of simple 
construction r containing no movable parts. The major components 





Fisfure 2.1. A typical solid propellant motor 



include the combustion chamberf which contains the propellant 
drain of particular deometry which in turn determines the surface 
burnind area. Combustion of the drain denerates hidh temperature 
dasesr which are ejected throudh an exhaust nozzler desidned to 
accelerate the dases to as hidh a velocity as possible* The 
throat is the point in the nozzle of least diameter^ and the head 
beind the front of the motor. Combustion is initiated by passind 
an electrical current throudh an idniter which contains a charde 
of black powder. 

The theoretical analysis of the rocket motor necessitates 
certain simplif icationsr that is» the assumption of an ideal 
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rocket4 The usefulness of this concept is indicated bw the fact 
that measured performance can be expected to be within 1 to 10 
per cent of calculated ideal values C23* 

An ideal rocket assumes the followina* 

(1> The propellant combustion is complete and nonvariant to that 
assumed b« the combustion eou3tion4 

(2> The workinsl substance obeus the perfect das law4 
<3) There is no friction4 

<4) The combustion and flow in the motor and nozzle is 3diabatic4 
<S) Steady-state conditions (unless stated otherwise) 4 

(6) Expansion of the workind fluid occurs in a uniform manner 
without shock or discontinuities4 

(7) Nozzle flow is one dimensional4 

<8> The das velocitw> pressure» and density is uniform across 
ana cross section normal to the nozzle 3xis4 
(9) Chemical ectuilibrium is established in the combustion chamber 
and does not shift in the nozzle4 

Ana additional assumptions will be stated as necessara4 

242 Propellant 
24241 Composition 

The important considerations for amateur rocket propellants 
are the availability of the constituentsr cost? safety of hand- 
lind castabilitar consistency of performancer and adeauate perfor- 
mance4 The propellant that is investidated here easily meets or 
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eKceeds these reeui rements* 

The propellant consists of sucrose (susfar) fuel with potas- 
sium nitrate as the oxidizer. The 0/F ratio is chosen on the basis 
of which ratio dives the greatest overall perfrmance for a diven 
propellant drain desidn. For most propellants includind this one» 
it is on the fuel rich side of stoichiometric. More accuratelur 
varuind the 0/F ratio affects both the characteristic exhaust 
velocity* CV* (essentially relates to the enerdy available)* the 
burnrate* r* and the Quantity of non-daseous exhaust products. 

These aspects are treated in closer detail elsewhere* 

2.2.2 Combustion 

The assumed combustion eouation is based on the reaction of 
potassium nitrate (KNOg) with oxyden ^02^ C33* and upon the decom- 
osition of sucrose ^^12^22^11 ^ heatind C43. The assumed 

combustion eouation is. 

a C^2^22°11 ^ »<N 03 “-> c CO + d CO^ + e HgO + f KgCOg + d Nj + 

h C + i Og 

where KgCOg represents potassium carbonate* where the coefficients 
a throudh i are dependant upon the O/F ratio with c*h or i often 
beind zero. 

This eauation assumes that secondary exhaust products (such 
as NO*K20» etc.) are formed in nedlidible Quantities. As well* the 
effects of dissociation are assumed nedlidible. This second assum- 
ption is denerally valid for combustion temperatures below 1700 K* 
particularly at hidher pressures C5D. It should be noted that 
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there are two non~aaseous products fortned» K_CO_ and carbonr from 

Z o 

which no expansion work can be derived. 

Knowledge of the combustion eouaiion allows the calculation 
of the adiabatic flame temperature (AFT)» the maximum possible 
combustion temperature. The AFT is determined from an enthalps 
balance! 



Z n f h“ + Ah = n h“ + Ah 

iLf Ji ^eLf 



< 2 . 2 ) 



where R and P refer to the reactants and productsf respectively f 
n. and n_ are the individual reactant and product molal numbers 

I ^ 

respective! w » Ahj » the enthalpy of formation per molei>Ah » the 
standard enthalpy at the specified temperaturer per moler and 
where! . 

f- 

Ah = / C dT + Ah 

\ P tr 

and Ah = enthalpy of transition/mole 

tr 

T as reference temperaturex typically 300 K 
1 

T « AFT 

K 

C as specific heat at constant pressure/mole 

p 



Appendix A contains a sample calculation for the case of 
0/F “ 65/35. 

Fidure <2.2> shows a comparison of theoretical to actual 
flame temperatures C73. 
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Figure 2*2« The relationship between calculated and actual 
temperature* 

The AFT varies with the 0/F ratio? particularly strondly 
in the rande close to stoichiometric? where the AFT is maximum? 
as shown in fidure (2«3)» 

The two points of sharp slope chande are at stoichiometric 
mixture (peak) where the 0/F ratio is 73»3/26»7 and the combus- 
tion eouation is 

^12^22°11 *^2 ^ ^2° 

(2.3) 

and the point where partial oxidation of C to CO is completed? 
with an 0/F ratio of 63.9/36.1 ? with the combustion eeuation. 
‘^12^22°11 ^ 9 CO + 11 HgO + 3 K^COg + 3 



(2.4) 
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Fiaure 2»3. Adiabatic Flame Temperature as a function of 0/F ratio. 

The 0/F ratio that is normally employed for this motor is 
65/35. This dives the followind combustion eouation. 

^12^22°11 KNOg— > 6.29 CO + 0.56 COj + 11 HjO + 

3.14 KgCOg + 3.14 N 2 <2. 5) 

The performance analysis of the rocket motor is based upon 
this combustion ratio* unless stated otherwise. 

Another parameter that is used extensively in ideal perfor- 



mance analysis is the averade molecular weidht of the daseous 
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exhaust products* This is calculated from the combustion 



products* 



and 


from 


the 


individual 


molecular weidhts of 


the 


M' = 


Jll. 


M' 


+ M' 


+ — M' + *** 


(2.6) 






i 


n^ J 


n^ k 





where ir J and k are the individual constituents* n is the mole 
number* and where t refers to the total number of moles* 

The ratio of specific heats* k » is another important 
Quantity in the analysis of compressible fluid flow* which is the 
type of flow encountered in rocket nozzles* where k is defined by 



P 

k = <2*7) 

where * for an ideal eas* k is a function of temperature only* 
The value of k can be determined from knowledge of the 
specific heats* Cp * of the individual exhaust eases* where 



k * (2*8) 

Cp - R' 

where R' is the universal das constant* and 




However* complications arise in the flow throudh a supersonic 
nozzle where the temperature of the eases drop appreciably (as 
will be shown later)* since Cp can be a strond function of temp- 
erature* It is possible to analyse such a flow utilizind a vari- 



able isentropic component C8I1* however* a sufficiently accurate 




11 



result can be obtained bu calculatind an average value of k for 
flow through the nozzle (see appendix B>» Testinsi has found that 
for a 15 “ conical nozzle with a small area ratio (6 to 8)» the 
variation in k does not have a pronounced effect C93« This is 
similiar to the nozzle tape used for the motor under consideration* 
The remainins combustion parameter to be determined is the 
burnraten r (also called the surface regression rate). This is 
normallH determined empiricallar and is a function of the propel- 
ant composition and certain conditions within the combustion 
chamber* These conditions include propellant initial temperature* 
chamber pressure* and the velocita of the Saseous combustion 
products over the surface of the solid (erosive burning)* It is 
necessary to combine a theoretical model toSether with empirical 
data in order to particularize the burnrate for a siven propellant. 
The usual model is to approximate the burnrate as a function of 
pressure* 

n 

r - a (2*10) 

where a and n are empirically determined constants* and is the 
combustion chamber pressure* The pressure exponent n» associated 
with the slope of the pressure-burnrate curve is almost independent 
of the propellant temperature* The coefficient 3 is a function of 
the initial propellant temperature* but not of pressure* From 
eouation (2*10) it can be seen that the burnrate is very sen- 
sitive to the exponent n* Hish values of n sive a rapid chanse of 
burnrate with pressure* This implies that even a small chande in 
chamber pressure produces substantial changes in the Quantity of 
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hot sl3ses produced* As n approaches unitur burnrate and chamber 
pressure became very sensitive to one another and disastrous rises 
in chamber pressure can occur in a feu milliseconds* Uhen the value 
of n is low and becomes closer to zerar burninsf can become unstable 
and may even extinguish itself* Most producton propellants have a 
pressure exponent randind from 0*3 to 0*6 C103* 

A correction is denerally made for erosive burnind where 
the increased burnind can be accounted for by an empirical 
correction of the form* 



r*r„<l+ku> <2*11) 

where k is an empirical constant and u is the das velocity. 
However » for this report this form of correction will not be 
applied* Insteadr the burnrate is calculated with the erosive 
term taken directly into account* 

2*2*3 Propellant Grain 

The drain is prepared by castind the molten propellant into 
the desired shape* Several common conf idurations are shown in 
fidure (2*4)* The main criterion for choosind the deometry is to 
achieve the desired thrust- time characteristics* Thrust is a 
function of the instantaneous burnind area which itself is de- 
pendant upon the initial drain conf iduration* The most common 
desidns are those which achieve prodressive» redressive» or 



neutral thrust- time curves^ as shown in fidure (2*5) 
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Fislure 2«4« Several common ^rain configurations* 




Fiaure 2*5* Internal burnins arain desians with their thrust-time 
characteristics « 
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The ^rain confiauraton of the motor beind considered is an 
internal- external burning tube» unrestricted at the ends» The 
theoretical thrust- time characteristic is slidhtly to moderately 
redressivei- deeendins upon the lensith/diameter ratio. 

The initial burning arear is siven bw. 

Aj^ = 7 t L (D^ + D. ) + 0.5 <Dg - dJ >] <2.12) 

where L is the drain lendth» and D. are the drain outside and 
inside diameters^ respectively. 

A slidhtly modified version of this drain is also used in 

the motor (fidure 2.6)» where A. is diven by. 

b 



A = 7T I (D^ + D,) + — -r3D„ <D + D^)l 

b I 4 L 2 



[t 



+ B, ) + - 2D, 



]| 



<2.i3) 



These two expressions can be modified by utilizind eeuation 
<2. 11) to obtain the instantaneous burnind arear A^>(t)« 

P * 

A,,(t) = 7T (L- 3 t + D, )+ 0.5 ( Dj, - D, ) I (2.14) 

and 



A (t) 
b 



7r|(L, 



n 

atP^ ) (D + D > + 
*» o 

1 







+ 





+ D ) 
i 




<2.15) 
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where t represents the time from initial burnins. It is important 
to recoSni^e that these expressions represent ideal burnins 
which assumes that simultaneous isnition of all surfaces occurs 
at the beainninsi of the burn* As well the effects of increased 
faurnina surfaces due to flaws in the arain such as bubble holes 
and other flawsr as well as the effects of erosive burnina are 
nealected» 

These expressions are useful in calculatina the chamber 
pressure and therefore thrusts as will be shown subseauently ♦ 
Propellant densitsu pp and arain density* Pg are two 
additional properties that will prove useful* Ideally these two 
are identical but as a result of voids in the arain its density 
is somewhat lower* 

The ideal propellant density is a function of the 0/F 
ratio* since the density of the two constituents are different 
(P^^o “ 2*11 a/cm* til]* '^sucrose" a/cm* C121 )♦ The ideal 

3 

propellant density can be expressed ast 
1 f o f f . ___ , 3 

— = — + — = i.888g/cm® <2,16) 

p9 P o pi 

where f and f. refe^r ^6“ tFe mass fraction of the oxidizer and 

O ff 

fuel* respectively* For the 65/35 0/F ratio the actual propellant 
density has been found to be about 5 percent lower than ideal* 
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CHAPTER 3J NOZZLE THEORY 

The snalasis of a rocket nozzle flow involves the studa of 
steadaf one-dimensional compressible flow of an ideal aas» The 
actual flow differs somewhat from this simplified model particu- 
larla in resard to the presence of solid or liauid particles in 
the flow stream* The necessara modifications will be dealt with 
in the next chapter* 

The analasis of compressible flow involves four eeuations of 
particular interest! continuita? momentumr enerSaJ and the 
eeuation of state* These eouations are applied to desiSn a nozzle 
with the objective of accelerating the combustion dases <and 
particles) to as hidh an exit velocita as possible . This is 
achieved ba desidnind the necessara nozzle profile with the 
condition that isentropic flow is to be aimed for * This 
necessitates minimizind frictional effectsf flow disturbances r 
and conditions which can lead to shock losses* As wellr heat 
losses would have to be minimized* 

3*1 Nozzle Flow 

In describind the state of a fluid at ana point in a flow 
field it is convenient to emploa the stadnation state as a 
reference state (the state characterized ba a condition of zero 
velocita )* Local isentropic stadnation properties are those 
properties that would be attained at ana point in a flow field 
if the fluid at that point were decelerated from local conditions 
to zero velocita followind a frictonless adiabaticr that is» 
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isentropic process. 

The energy eouation for an adiabatic flow between points 
X and y in which the decrease in enthalpy is eouai to the increase 
in kinetic enersta is Siven asJ 

1 2 2 

h-h - — (V - V ):= Cp<T-T ) <3.U 

X a 2 a X x a 



where h? v» and T are the enthalpy^ velocitar and temperatures 
respectively. 

The stagnation temperature is defined from the enerda 
eeuation as? 

2 

= T + V <3. 2) 

2 Cp 

For an isentropic flow processs the following relations for 
stagnation conditions hold? 




<3. 3) 



The local acoustic velocity of an ideal das is defined as! 



'a = k R T (3.4) 
where R is the das constant. The Mach number is defined as the 
ratio of flow velocity to the local acoustic velocity: 

V 



M 



a 



<3. 5) 
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From eauations (3.2) r (3.4) » and (3.5) the total temperature 
Mach number relationship can be written* 



T k-1 2 

_o = 1 + M (3.6) 

T 2 



It can be shown from the first and second laws of thermo- 
dynamics C133 that for an ideal das underdoinsl an isentrppic 
process assuming constant specific heatsr that* 

P 

-~j- = constant (3.7) 

P 

From this result and from the eouation of statef P = pRT» 
the stasination pressured density - Mach number relationship 
can be expressed as* 




(3.8) 



(3.9) 



The use of eeuations <3.6)r (3.8>» and (3.9) allow each 
property (TrPf P > to be determined in a flow field if the Mach 
number and the stagnation properties are known. From the enerdy 
eouation for an adiabatic flow (3.1) the stadnation enthalpy is 
defined. 

2 



h = 

o 



h + 



(3.10) 
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Phusicallur the stagnation enthalpy is the enthalpy that 
Mould be reached if the fluid were decelerated adiabatically to 
zero velocity. We note that the stastnation enthalpy is constant 
throughout an adiabatic flow field* Since the above stastnation 
properties (P^ t and T^) are related to the stagnation enthalpy 
by the specific heats and by the eouation of states it is 
apparent that each of these stagnation properties is constant 
throusfhout the adiabatic flowfield* 

For steady? one-dimensional flow? the continuity eauation 
can be written 

* t * 

pAv = constant = P A v (3.11) 



where A is the passage area? v is the velocity of the flow? and 
a starred (*) variable indicates critical conditions? or where 
M is unity. 

From eouations <3.4>? (3.6)? (3.9) and (3.11)? it is 

* 

possible to express the area ratio? A/A? in terms of the Hach 



number 



A 



1 

H 



1 + 



1 + 



k-1 

2 

k-1 

2 



a{k-i) 



M 



(3.12) 



Figure (3.1) clearly shows that a conversins-diverains 
passade with a section of minimum area is reouired to accelerate 
a flow from subsonic to supersonic speed. The critical point 
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where M becomes unitu is seen to exist at the throat (point of 
minimum area) of the nozzle» 

The variation of the properties durinsl flow through the 
nozzle is illustrated in figure <3*2). From eauation (3*1 >i' 

the nozzle exit velocitaf can be found b«* 

o 



V 

e 




(3.13) 



Figure 3*1. 




Variation of A/A with Mach number in isentropic 
flow for k = 1.4. 



where h and v are the enthalpa and velocity at ana point 

X X 

within the nozzle. This relation applies to both ideal and 
nonideal rocket units. This eeustion can be rewritten with the 
aid of eouations (2.8) and (3.3) 






Fia. 3.2 . Msristion of density Pr pressure Pf snd 
temperature T» throudh the rocket nozzle. 

for flow between the chamberr where stadnation conditions are 
assumed to exist> and the nozzle exit. 

From this eauation it is seen that the maximum exit velocity 
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occurs 3t sn infinite pressure ratio Fj^ /P® » or when exhaustind 
into 3 vacuum* 

The ratio between the throat and anu downstream area at 
which pressure prevails can be expressed as a function of the 
pressure ratio and k as follows* usins eouations <3«3>» (3»9) *( 3 * 4)7 (5,4) 
and (3*14) (notinS that at the throat M is units )♦ 



* 

A 

A 

X 




<3*15a> 




This eouation is important in that it allows the exit area* 

A_^ to be determined such that the exit pressure is ectual to the 

e 

ambient pressure P. (tspicalls 1 atm«)« This is known as the 
desiSn condition where it will later be shown that for such 

conditions maximum thrust is achieved. For this desian* A / A* is 

e 

known as the optimum expansion ratio. 



3.2 Rocket Performance Parameters 

This section deals with the various performance parameters 
that are used to determine and compare the performance of solid 
propellant rocket motors. As well* modifications to the simpli™ 
fied model are considered to correct for real or actual perfor- 
mance. As well the effects of two phase flow are considered. 
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3.2«1 Thrust 

The thrust F» of m rocket motor can be shown to be siven 
bu C143: 



F=/pdA =mv +(P-P)A (3.16) 

e e 3 e 

where the first expression represents the intei^rsi of the pressure 
forces 3ctina on the chamber and nozzle projected on a plane normal 
to the nozzle axis (fidure 3»3>f m is the mass flowrate of the 

exhaust products and v^is the exit velocity* The second term of 

e 

the second expression is called the pressure thrust and is eoual 
to zero for a nozzle with an optimum expansion ratio* 

From continuityf eeuation (3*16) can be rewritten! 



» « 

F=PvA v+(P-P)A 

e e a e 



(3.17) 



and modified usin^ (3*14)r (3*9)» and (3*4) to yield* 




(P - P )A 

0 B 0 

(3.18) 



where the eouation assumes k is constant throughout the expansion 



process* This eouation shows that thrust is proportional to* 

* 

i> throat area? A 

ii) chamber presssurex P^ 

iii) pressure ratio across the nozzler 

iv) specific heat ratiox k 

v) pressure thrust 
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Fi£f» 3*3 « Pressure balance on chamber and nozzle wallsr 

internal aas pressure is highest inside chamber 
and decreases steadils in nozzler while external 
atmospheric pressure is uniform* 



3*2*2 Thrust Coefficient 

The thrust coefficient r is defined as the thrust divided 
by the chamber pressure and throat area* 



C, == ----- (3*19) 

Po A 

The thrust coefficient determines the amplif ication of the 
thrust due to the sas expansion in the nozzle as compared to that 
would be exerted if the chamber pressure acted over the throat 
area only* From eauation (3«18)« 




(3.20) 



For any fixed pressure ratio P® /P© the maximum C^ can be 
found by takind the derivitive dC^ =» 0* Thereforef the 

d(A /A* > 

0 
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maxintum Cj will occur when P « P » or when the nozzle is designed 

e 3 

for optimum expansion (desidn condition)* 

Eoustion (3*19) is useful for comparing the measured » 
to the theoretical value as determined from eouation (3*20)* 



3.2*3 Characteristic Exhaust Velocity 

The characteristic exhaust velocityr CV» is defined as 
CV c /Cf where c is the effective exhaust velocity* 



= V + ( p - p ) — z 

• * 

m e earn 



(3*21) 



The CV can be expressed as a function of the sas properties 
in the combustion chamber usind eauations (3.21) r (3*19) » (3*1 )»(?.9) 
and (3*11)1 

(3*22) 




CV is usually used as a fisfure of merit of a propellant 
combination and combustion chamber desistnr and is essentially 
independent of nozzle characteristics* This makes CV useful as 
a comparison for different propellants* 

T r recallins that 



It is interestind to note that CV 



H' 



R = R'/ M' * A hidh value of CV is desirable?, however a hidh value 

of T midht not be* Hidh combustion temperatures may cause ex- 
o 

cessive heatind of the parts of the rocket motor that are exposed 
directly to hot dases* Nozzle erosion can be severe with conven- 
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tlonal materials particularly when burn durations exceed several 
seconds# For this reason it may be more desirable to increase 
CV by reducina M'' rather than increasind « This can often be 
accomplished by usind a propellant that is on the fuel rich side 
of stoichiometric* 

The CV can also be considered an expression of the impetus 
(n R T) of the propellant# This is useful for determinind a 
propellant in the laboratory by means of a closed bomb 
(constant volume) measurement ClSl# 



3#2#4 Impulse 

The impulse (or total impulse) rlr is the intedral of the 



thrust over the operatind durationi» tJ 



I = 




( 3 4 23 ) 



or the area under the thrust- time curve (fidure 3#4)# 




Fidure 3#4» Area under thrust-time curve represents the 
total impulse of the motor# 
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The operstina duration is approximately eeual to the burn 

timex t r far motors with londer burn times (t. > several seconds) 
b ^ 

For short burntimes such as those associated with small motors 

the duration of thrust would be considered the burntime plus the 
time duration for the residual dases to exhaust from the combus- 
tion chamber after burnout* The total impulse would then be diven 
by 

2 

I = C, P- A* T <3*24> 

k+1 

where r is a time constant aiven by 

1 A*Po 
r yC7 o CV 

where n is the bufhrate exponent ahd“V is -^ne chamber volume* 

o 

The specific impulse is one of the most important perfor- 

• V 

mance parameters used in rocket research* It is defined as the 
thrust that can be obtained from an eouivalent rocket which has 
a propellant weiaht flow of unity. It is aiven by: 

F c 

I = --- = (3*25) 

sp w a 

With solid rocket motors it is difficult to measure the 
weiaht flowrate* so the averade specific impulse is usually 
employed * where I = 1/ t* 

SP 

The ideal specific impulse can be calculated for a aiven 
propellant and motor combination from eeuations <3*2S> arjd (3.14): 



» 
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<3»26) 

For comparison purposes is usu3ll\:t taken as 1000 psia 
and P^ as lA ,7 psiar b« convention* A simplified semi-empi rical 
method for determining specific impulse was developed b« Free 
and Sawuer C171 with a claimed accurac^:^ of 3 to S percent* 
Recalling that CO is also a measure of the propellant 
enerdsf the actual specific impulse can be determined bu utili- 
zing the ballistic bomb method as previously mentionedr where 
I is related to CM byj 

SP 

CM C, 



I = (3,27) 

SP d 



however^ reouirind the use of the ideal C^ unless the actual value 
is known* 

3*2*5 Chamber Pressure 

The combustion chamber of a solid propellant rocket is 
essentially a hidh pressure tank containind the entire solid mass 
of propellant* Combustion proceeds from the surface of the drain 
where the rate of das deneration is eoual to the rate of consump- 
tion of solid material r in the ideal case where only daseous pro- 
ducts result <no solids or liouids)* 

m = P. A, r (3*28) 

d 

where m is the das deneration rate? A. is the area of the burnind 

d ** 

surface» and r is the surface redression rate (burnrate)* The rate 
at which das is stored within the combustion chamber is diven by* 
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dirij d 

- ( Vp ) (3»29) 

dt dt 

where is the instantaneous Sas density and is the instantaneous 
chamber volume> The use of these two eouations as well as the 
expression for the rate at which Sas flows throusJh the nozzle and 
the expression for burnrate» leads to the fallowing expression 
for chamber pressure IZ183* 




For constant and therefore constant thrust it is clear 
that the burning area must remain constant» such as with a neutral 
burnina drain. The instantaneous chamber pressure can be diven by 
the same expression* in which is the instantaneous burnind 
area. The variation of A with time depends upon the burnind rate 
and the initial deometry or the propellant* drain* 

For short burn time motors such as the one under consid- 
eration (less than one second burntime) the duration under which 
steady-state conditions exist* where eeuation (3.30) applies* 
accounts for perhaps one- third to one- half of the total thrust 
time. Prior to idnition* the combustion chamber is filled with 
cool air of ambient pressure. The start-up period where the 
pressure builds up to the steady-state value is a complex 
problem that has been studied by Von Karman and Malina [!191* 

They have developed an expression for the pressure as a function 
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of time during st3rt-up» ^iven by* 



P = pQ ( 1 “ (3,31) 

where P is the instantaneous pressure at time t and where r is 
previously defined* 

Immediately after burnout the combustion chamber remains 
filled with hiah pressure aas. The expression for chamber pressure 
is siven as a function of time by C201* 



P 





(3.32) 



3.2*6 Corrections for Two Phase Flow 

The previous analysis of nozzle flow and performance para- 
meters considered the ideal case where the workina fluid is a pure 
ass* This analysis is fine for the case where there is no solid 
or liouid particles in the exhaust (such as for liouid propellants) 
but must be modified for the case where such particles are present* 
The presence of particles is detrimental to the overall rocket 
performance» with the extent dependant upon several factors* 

One of these factors is motor size* It has been found by 

4 

Gilbertr Allport» end Dunlop C21II that for larae motors < >10 

pounds thrust) the overall specific impulse losses are low)> but 

2 

for small motors ( «10 pounds thrust) the losses can be siani- 
f leant* The motor under consideration is of this maanitude thrust 



makina corrections for two phase flow necessary* 
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The losses in performance are due to two factors* 

i) velocity lasl» where the particles are not accelerated 
through the nozzle to the same velocity as the sasesJ and 

ii) thermal las» where the particles are not in thermal 
eouilibrium with the surrounding sases* 

It has been found that the latter has a small effect? usually 
allowinaf the assumption that the particles are in thermal 
eeuilibrium C221* 

The starting point in the consideration of the effects of 
solid (liauid) matter formation is with the eouation relating to 
the Sas Generation rate (eouation 3«28)» This eouation becomes* 



m = ( 1 ~ X ) r (3*33) 

sf " 

where X is the mass fraction of solid (liouid) in the exhaust* 

As a conseouence? this would chanae the eouation for steady-state 
chamber pressure toJ 

1 

1-n 

(3*34) 



where the chamber pressure is reduced by a factor of ( 1 - X > * 

The expression for specific impulse would chanGe to account 
for the difference between particle and Gas velocity* 




I 

SP 




(3.35) 



where i denotes a particular constituent* If all the particles 
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ere assumed to have the same velocity this can be expressed as* 



1 <m V + m„ v^ ) 
s s g g 



<3*36) 



+ "-3 



sp s 

To accurately describe the effects of particles on a 
rocket's performance rectuires knowledge of heat transfer proper- 
ties and dras processes of the particles in order to express v^ 
and Tj in terms of the overall flow properties* It is possiblej 
howeverf to derive the limiting cases to determine the least and 
most detrimental effects* 

In the first model r the particles are assumed to be very 
small* where T « T_and « v_ * This would dive the followind 
expression for exhaust velocity (see appendix c for derivation)* 




(3*37) 



where Cj is the averade specific heat of the solid (liouid) matter* 
Note that this eouation reduces to the familiar form (3*14) when 
X = 0 (no particles)* 

The expression for thrust would be modified to this form 
(see appendix D for derivation)* 



A P -- 




+ (P - P ) A 
e a e 



(3.38) 
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Note 3^3in that the eouation reduces to the familiar form <3«18) 
when X = 0» 

The expression for the thrust coefficient is derived in the 
same manner^ t:<ieldins the expression? 




It is interestinst to note that the thrust coefficient in- 
creases with the presence of particles compared to the slas only 
value* This eouation implies a 75 percent increase in the thrust 
coefficient for a typical propellant where X = 0.45* 

The second limitins case models the particles such that the 
particle velocity and temperature remain essentially constant* 

This implies a combination of relatively larsJe physical sizer 
slow heat transfer^ and low dras* With such a model » the expres- 
sions for the performance parameters remain identical to the 
original (das only) expressionsr i*e. the das flow is not affected 
by the presence of particles* Howeverr the specific impulse is 

reduced* From eouation (3*36) where m v « m v * 

s s d d 



1 

I - ( 1 - X > V <3, 40) 

sp d d 

Thereforer the specific impulse is reduced by a factor of 1-X* 

Research conducted into two phase flow also indicates that 
the nozzle type has a bearind on the effect of particles in the 
exhaust stream T231* As well rit has been found that particle 
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size (and distribution) is independant of motor size» and is 
therefore larSela a property of the particular propellant C243. 
Another detrimental effect of liouid particles is the tendency 
for the material to build up on the nozzle surfaces? reducing the 
effective nozzle area* Additional losses are introduced ba the 
resulting rousfh surface which increases skin friction C25J* 

For the motor beinS considered? the losses due to particle 
presence will be accounted for in the enersa conversion factor? 
discussed below* 

3*2*7 Corrections for Real Nozzles 

The precedina analysis considers ideal rockets? which of 
course do not exist* The ideal case represents the maximum per- 
formance that can be attained? with the actual performance beind 
reduced ba a number of factors* These factors are accounted for 
in real nozzles ba usinS various correction factors* 

Conical nozzles reauire a factor to correct for the non- 
axial component of the exit aas velocity as a result of the 
divergence ansle? defined as 2 a (where a is the half ansle)* 

This factor N ? is Siven ba* 

1 

X = ( 1 + cos a > (3*41) 

2 

As well? the flow in an actual nozzle differs from ideal 
because of frictional effects? heat transfer? imperfect sfases? 
nonaxial flow? nonuniformita of the working fluid and flow 
distribution? and the effects of particles* The desiree of de- 
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parture is indicated ba the eneraa conversion efficienca of the 
nozzle. This is defined as the ratio of the kinetic enersla per 
unit of flow of the Jet leavinsi the nozzlei- to the kinetic eneraa 
per unit flow of a hapothetical ideal Jet leaving an ideal nozzle 
that is supplied with the same working fluid at the same initial 
state and velocita and expands to the same exit pressure as the 
real nozzle. This is expressed ast 



e 




(3.42) 



V 



where the subscpipts i and a refer to the ideal and actual states 
and e denotes the enerau conversion efficiency. 

The velocity correction facton>r, » is defined as the sauare 
root of e. For most production motors» the value lies between 0.85 
and 0.98. This factor is also approximately the ratio of the 
actual to ideal specific impulse. 

The discharge correction factor r is defined as the ratio 
of the mass flow rate in a real nozzle to that of an ideal nfczzle 
that expands an identical workind fluid from the same initial 
conditions to the same exit pressure? 

m me 



" m. F 

I 

The value of the discharge factor is often 
unityf because the actual flow can be larder than 
for the followind reasons? 



<3. 43) 

larder than 
theoretical 



i) the molecular weidht of the dases increases slidhtly 
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when flowing through s nozzler thereby chanainsi their density* 

ii> some heat is transferred to the nozzle wallsr lowering 
the das temperature » increasinst its density* 

iii) the specific heat and other das properties chande in 
an actual nozzle in such a way as to slidhtly increase the value 
of the discharde coefficient* 

iv) incomplete combustion increases the density of the 
exhaust dases* 

These correction factors result in a thrust lower than for 
the ideal case* 



F =5 
a 



f X F 

d ; 



F 

a 



f 



V 







< 3 * 44 ) 



( 3 * 45 ) 
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CHAPTER Ai EXPERIMENTAL TECHNIQUE 
4*1 Motor 

The motor under consideration is illustrated in the fiaure 
below (fiaure 4*1) » with a complete description aiven in appendix 

< E ), 



Z5 cm. 




Fiaure 4.1* Actual motor used in testina* 

The standard nozzle has a conical profile with a 12 dearee 
diveraence half snale* It was fabricated from mild steel bar stock 
turned and bored to the desired dimensions on a standard lathe* 

The nozzle contour is rounded at the throat to avoid sharp dis- 
continuities which could lead to shock losses* As well* the entire 
nozzle inside surface is polished to reduce friction losses* 

The nozzle is retained to the combustion chamber bw six hiah 



strenath alien screws 
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The head of the motor is readily removable to allow for 
easy loading of the strain* Constructed as well from mild steel f 
the head is retained by a sinstle safety bolt (3/16 X 1 1/2 in.? 
grade 6) designed to shear if the chamber pressure should attain 
an unsafe level (it MPa*)* The combustion gases are sealed from 
escaping around the head by five asbestoes-f ibre composite gas- 
kets* The nozzle is sealed against the combustion chamber by 
utilizing a tool which rolls a circular die under high pressure 
around the perimeter of the chamber outside wall against the 
nozzle section forward of the retention groove* This effectively 
seals the chamber wall against the nozzle to reduce gas leakage* 

4.2 Cubic Nozzle Profile 

Although the exact nozzle profile is not critical for 
good performance as the flow occurs in a region of favourable 
pressure gradient* optimum performance is achieved by designing 
the nozzle by the method of characteristics * A thorough dis- 
cussion of this technioue is presented in references C263 and 
C273* This method is lengthy and complex and will not be consid- 
ered here* However* an approximation to this design is that of a 
cubic profile where the nozzle contour follows the curve of a 
cubic eouation C28D* 

A nozzle was constructed with the aid of a computer numer- 
ically controlled lathe. The criterion for design was based on the 
same throat* entrance and exit areas as the conical nozzle* The 
convergent section was based on a maximum half angle of 30 degrees* 
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This nozzle was constructed of cold rolled mild steel 
bsr stockr end was highly polished on the flow surfaces> The 
diverslence exit ansile was machined to zero dearees» to eliminate 
radial flow losses <non-~truncated) ♦ 

It should be noted that while the critical dimensions 
of the cubic nozzle are coincident with the conical nozzleo the 
overall lendth is about 17 percent shorter. 

4.3 Propellant Preparation 

As stated previously? the standard propellant com- 
bination is 65 % oxidizer <KN0g) and 35% fuel (sucrose). The 
propellant (srain) preparation techniaue for any 0/F ratio is 
identical? however. 

In order to ensure complete mixinst and intimate contact 
between the oxidizer and fuel crystals? each constituent must 
be sfround to a small particle size. Commercial oxidizer crystals 
are around according to the followina standard ranaest 

—6 

Coarse 400 to 600 M ( 1 *< = 10 m) 

Medium 50 to 200 M 

Fine 5 to 15 

Ultrafine submicron to 5 M 
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The commercially obtained potassium nitrate oxidizer has 
an average particle size of 250 The oxidizer is subseauently 
dround down to an averasie particle size of 100 1* ♦ The sucrose is 
obtained in powder form with typical particle size of approxi- 
mately 10 M » 

It has been found that the particle size of the oxidizer 
has an important influence on both the burnind rate and more 
importantly on the impulse delivered* Therefore? careful attempt 
was made in the experimentation to obtain consistency in oxidizer 
particle size* 

After accurately weisihind the desired quantities of both 
constituents? they are placed in an electric rotatinsf drum mixer 
where mixind occurs for twenty hours* 

The drain castind procedure involves heatind the mixture 
in a container maintained at a temperature of 190 to 200 dearees 
celcius by immersion in an electrically heated oil bath* The 
sucrose? <m«p* 182 dedrees C*> combines with the non-meltind KNOg 
(m*p» 441 C*)? to form a slurry* Decomposition of the sucrose 
initiates as well at this temperature beind both a function of 
of temperature and time* Therefore? the temperature has to be 
maintained within this rande? and for the minimum time duration 
necessary for the entire mixture to melt* 

This decomposition (caramel ization) ? if kept in check? has been 
found to have no detrimental effect on the propellant performance* 
Investidation of the effects of low dedree caramelization of 
sucrose indicates that there is no sidnificant mass chande? which 
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is irtiportent since the occurance of such would alter the effective 
0/F ratio. As well» it has been found that less than 1 percent 
moisture is absorbed b« the finely divided suSarx which is asain 
an important consideration for any moisture initially present 
would be diven up upon heatins* 

Once the slurry has become sufficiently fluidx it is 
cast into the desired shape* This shape (hallow cylindrical) 
is achieved by pourinsi the slurry into a lubricated cylindrical 
mould* of slishtly smaller diameter than the combustion chamber. 

A lubricated bore rod is inserted down the central axis. After 
allowind to cool and harden for approximately forty-five minutes* 
the srain is removed* trimmed to size* weisfhed and measured* The 
strain is then stored in a sealed container. This is to prevent ex- 
posure to the open air for the propellant in its final form is 
very hydroscopic. 

4.4 Oaryind 0/F Ratio 

A number of drains were prepared to determine the effect 
of varyind 0/F ratios upon the motors impulse. As well* several 
small cylinders (1.5 x 10 cm) of varyind 0/F ratios were prepared 
for burnrate testind. 

4.5 Burnrate testind 

In order to Quantify burnrates under conditions of one 
atmosphere and at room temperature* several samples of different 
0/F ratios were prepared. The samples were cast into cylinders 
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3S indicated earlier* These cylinders were mounted uerisht by 
adherinsi one end to a base* A specified lendth was marked off» 
typically 7 cm* The top end was then isinited by a flat surfaced 
tool» heated prior in a flame* As the sample burned downx the 
elapsed time was measured between marks* This type of burnrate 
testins dives results for conditions of constant temperature 
propellant and non-erosive burnind* 

More complete burnrate testind would involve determinind 
the dependance of burnrate on pressure* In order to conduct such 
testind » it would be necessary to construct a suitable vessel 
that would allow burnind to occur at constant pressure* Such an 
investidat ion has not yet been conductedr however* Research of 
this type has indicated that the burnrate measured by this tech- 
nioue is not in complete adreement with the actual burnrate 
that occurs in the actual motor» beind about 7 percent hidher 
C29]* No eKPlanation could be found for this discrepancy* 

4*6 Motor Static Testind 

Static testind of the motor (fidure 4*3) to obtain the 
thrust-time characterics is conducted on a specially built test 
rid <fidure 4*4)* A remotely located Heath H-8 microcomputer is 
used for data acauisitionr with the measured sidnal sent to it 
via a 50 metre shielded cable* The acauisition system is potent- 
ially capable of handlind up to four channels of data input (e.d* 
thrustr chamber pressuref temperatures and wall stress)* For this 
series of testsr only thrust measurements were attemptedp however. 
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Fisiure 4.3« Static firing of the rocket motor mounted 
in test ria» 



4«6«1 Test Rid Construction 

The frame of the rid is constructed of heavy steel I-beam 
to eliminate undesired deflections under high thrust. Although 
the maximum thrusts endured by the motor under consideration 
uere 1.1 kN (250 lbs.)» the rig can handle much larger values of 
thrust . 

The motor is mounted vertical ly» nozzle upwards in a 
holder designed to allow for the small vertical motion encountered 



by the motor during firing 
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1 * 


Deflection bar 


7. 


Electrical idniter 


2. 


Lower stoe 


8 » 


Deflection bar tensioner 


3. 


Upper stop 


9* 


Frame 


4* 


Strain transducer 


10. 


Hydraulic damper 


5. 

6» 


Motor mount 
Rocket motor 


11. 


Damper adjust valve 



Figure 4*4* Rocket motor static testing risl« 
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The frsme supports a double cantilever deflection bar 
assemblu< It is adainst this bar that the motor acts durind firind 
resultind in the bar deflectind* The thrust is therefore trans- 
duced into a lineal deflection* In order to obtain hidh resolution 
of the thrust transduction* it is desirable to have a relatively 
larde defection of the bar* This maximum deflection is determined 
b« the transducer unit* havind been determined to be 0.635 cm. 

(1/4 in.). By knowind the maximum expected thrust of motor it is 
possible to choose the bar-cantilever assembly to achieve this 
(for any motor). 

The deflection y for a double cantilever assembly is diven 



by C303. 



3 

F L 

y = 

192 E I 



(4.3) 



where F is the force (thrust)* E is the Yound's modulus of the 
bar material* I is the moment of inertia of the bar* and L is the 
bar lendth. For a rectandular bar of thickness d and width b* this 
can be rewritten ast 

3 

F L 



y = 

3 

16 E b d 



(4.4) 



The deflection is therefore proportional to L/D cubed 

and inversely proportional to the width. Therefore* the choice 

of the bar dives dreat flexibility in achievind the maximum 

deflection* y . 

max 



A lower stop is placed directly below the centre of the 




bar to prevent damage to the transducer unit in case of over- 
thrustindu or a blowout of the motor head (see fidure 4»5)« 

The thrust induced displacement of the transducer bar 
results in an eaual displacement of the end of the transducer 
unitr located directly under the motor mount* This transducer 
unit uses a strain daude bridde fastened to its fixed end* This 
strain daude bridde forms the fundamental component of the data 
aouisition systenif illustrated in fidure (4*6)* The amplifier 




Fidure 4*5* Deflection bar assembly of static test rid* 




Fidure 4*6* Data aceuisition system used to measure 
motor thrust-time characterics * 
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circuit is reproduced in Appendix < F )♦ The A/D converter is an 
8“bit» four channel unit capable of eaniplinid 581 points per second 
(samples every 1«72 milliseconds)* These points are stored in the 
computer memory for Post-processinS ♦ 

Due to the sprina-mass nature of the deflection bar assemblyr 
it was found necessary to provide adeouate dampinsi to reduce 
oscillations of the deflection bar* A variable hydraulic damper 
was designed and built to allow a larde decree of flexibility of 
damping* The damper is detailed in fislure (4*7). The operation of 




SCREW VALVE 



Fiaure •4.7* Hydraulic damper used with test ria. 

the device is straiahtforward* A connectina rod with a piston at 
one end is attached to the deflection bar. The reversed vertical 
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motian of the rod alternately forces fluid from and back into the 
fluid chamber through a small Jet* The effective area of the Jet 
is controlled ba an adjustable needle valve» the adjustment of 
which doverns the dedree of dampina. This allows for a wide rande 
of dameind from essentially nil <needle turned out) to complete 
dampind (needle turned in)« 

The dampind fluid is standard hydraulic fluids and the 
system is bled after fillind to remove all traces of air. 

4 . 4.2 Calibration 

Calibration of the static test rid is performed with the 
aid of an e«tention arm with a series of weiahts hund at its end. 
The arm effectively amplifies the force which is applied at the 
deflection bar by a factor determined by the arm lendth (as is 
illustrated in fidure 4.8). 




Fidure 4.B. Extension arm setup for calibration of the 
test rid. 
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The force applied at the deflection bar F » is sliven byJ 

F =* — — - + U (4«5> 

L. ^ 2 

where U is the total weisht huna at the bar endf is the arm 

a 

weiaht« Since U^, r y Lj are constant for a aiven armi this 
can be written asi 

F = C,W + Cj <4«A) 

where = L^/ L| and /2 ♦ 

Calibration is carried out by hanaina a series of weiahts and 
recordina the output voltaae» The calibration data is presented 
in the form of voltaae V» as a function of forcer havina fitted 
the points throuah a second order curve <throuah the oriain) 
to yield a calibration curve of the followina form* 

2 

F = CgO + C^V (A. 7) 



where C* and C. are constants 

3 4 
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CHAPTER 5t RESULTS (Theoretical and Actual) 

5.1«1 Nozzle Testing 

The thrust-time curves for the conical and cubic nozzle 
are presented in figure (5«1) and (5»2)i respectively* The time 
period between successive points is 8*6 milliseconds <i*e. every 
5th point is printed) « For comparison purposes^ both results are 
presented on a sinstle curve illustrated in fisure (5.3) (notet 
different time base). 

A summary of the performance for both of these nozzles 
is presented in Table 5.1. 





Conical 


Cubic 


Total impulse 


(N-s) 


288 


281 


Specific impulse 


( s ) 


130.5 


127.3 


Max. thrust 


( N ) 


1155 


1075 


Thrust duration 


( s ) 


0.38 


0.40 


Propellant mass 


(kd ) 


0*225 


0.226 



Table 5.1 Comparison of conical and cubic nozzle motor 
tests. 



Both tests were conducted under similiar environmental 
conditions? although on separate dates. 

It should be noted that neither nozzle suffered any 
erosion or other detrimental effects after firind. This has 
been found to be the case even after multiple firings ( > 10) 

making the motors indefinitely reusable. 
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FiiSure 5.1. Actual thrust-time curve for the conical 
nozzle motor test. 





THRUST 
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Figure 5*2 



Actual thrust-time curve for the cubic nozzle 
motor test* 



74.6750031 

27.7819996 
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Fiafure 5»3. CotiiF>3rison of the conical and cubic nozzle 
motor thrust-time curves* 

5*1.2 Motor Testinsf of Varuind 0/F Ratio Grains* 

The results of this series of tests are presented in 
Table 5*2* Only the values for the average specific impulse are 
presented because it is felt that this would provide the most 
useful base for comparison* It should be noted that a non- 
standard conical nozzle was used for this series of testsf result- 
ing in lower overall values of specific impulse as compared to 
the standard conical nozzle as described earlier* Thereforer the 
results should be used for comparison purposes only* 



The impulse for the 75/25 0/F ratio was not measurable due 
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to poor burning resulting from the high oxidizer percentage* A 

large Quantity of liauid matter (potassium carbonate) was ejected 
from the nozzle» and burning continued for about half a minute* 



0/F Ratio 


a 

1 


50/50 


49*9 


55/45 


93*6 


60/40 


100 


65/35 


106 


70/30 


106 


75/25 





Table 5*2 Comparison of varying 0/F ratio motor tests* 

The 50/50 0/F ratio grain also resulted in poor firing. 

Large volumes of carbon were ejected with much remaining in the 
combustion chamber after firing in the form of a single porous 
mass* 

Each of these test firings were recorded on 8mm colour 
movie film for post firing examination of the test* The results 
showed that the exhaust flame varied in size and colour with 
differing 0/F ratios* For the more fuel rich ratios^ the flame 
was smaller and more orange in colour* For the more oxidizer 
rich ratios» the flame increased in size with a more purple 
coloured hue» a result of the potassium compound in the exhaust* 
As wellf it was found that the nozzle throat became hotter with 
increasing oxidizer percentager with the exception of the 75/25 
0/F ratio* The throat colour ranged from a dull red (50/50 0/F 
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ratio) to a bright uellow (70/30 0/F ratio)* This result is 
consistent with the theoretical prediction for combustion temp- 
erature as a function of 0/F ratio as discussed earlier* 

A comparison of the propellant properties T^i» H'f CO 
and X for varwins 0/F ratios as predicted bu theory is shown in 
f isure (5*4)* 




40 44 48 5T~ 56 60 64 68 72 76 ^ 



PERCENT OXIDIZER 



Fiaure 5*4* Comparison of theoretical properties! 

(1) Combustion temperature» 

(2) Averase molecular weiaht of exhaust 
aasesr 

(3) Characteristic velocityf CM 

(4) Particle mass fraction» X* 





5»i*3 Burnrate Testinsi 



The results of burnrate testing for various 0/F ratios 
are presented in fiiSure <5*5)» The conditions of testind were 
atmospheric pressure and room temperature* The burnrate results 
for the 75/25 0/F ratio are not shown due to the erratic burning 

encountered* 

As mentioned earlierr burnrate testing has not uet been 
conducted for conditions of elevated pressures to determine the 
burnrate-pressure relationship. Howeverr it is possible to est- 




Fislure 5*5* Comparison of atmospheric burnrate testing 
for varyinsl 0/F ratios* 

imate the pressure dependence in order to determine an approximate 
value for both the coefficients a and n (from eouation 2*10) » 
usins the results of the motor testing. This method involves 
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estiniatinM the time reouired for the drain web (wall thickness) 
to burn throudh under motor operatind conditions* The validity 
of this method lies in the assumption that essentially all the 
surface redression occurs at hidh pressure* This pressuref taken 
as P * occurs durind the steady-*state operation of the motor» 

O 

allowind an estimation of this time period from the thrust-time 
curve (see fidure 5*4)* 




Fidure 5*6* Estimation of web burnthroudh time 
from actual thrust-time curve* 



Since burnind prodresses from both the inside and out- 
side surfaces of the drainn the web ewpension rate is effectively 
doubled* The burnrate at a pressure of can then be expressed 



as « 




2 t^ 



P--Po 



n 



(5*1) 







59 



where end R. are the j^rain outside and inside radii« 

Combining this eeuation with the results of the atmos- 
pheric burnrate testinst 

r, = a <5»2) 

'P = P«.m 

and the steadu-state pressure eeuation (3*34) » results in a set 
of three simultaneous nonlinear eauations in an ni> and P * 

At a pressure of one atmospheren the burnrate is taken as 
r = 0*335 cm/s (65/35 0/F ratio)* Utilizing the thrust-time 
curve as shown in fisure (5*1) » the burnrate at can be approx- 
imated at 2*22 cm/s* Eouation (3*34) is applied with the 
following values shown below* 

3 

pp = 1*83 s/cm 
2 

A « 280 cm 

r 2 

A* » 0*6825 cm 
k « 1.283 
T « 1620 K 

The value for A, stated above is 10 percent sreater than 

b 

the value calculated by eauation (2*13). This is an estimation 

of the actual surface area with surface flaws taken into account* 

Solutions of eouations (5.1) » (5*2) and (3*34) with the above 

values dives the followinst results for a? n and P^ • 

a - 0*85 

n = 0*407 

= 10*5 MPa* (1530 psia) 

When the coefficients above are inserted into the ex- 
pression for burnrate (eauation 2*10)r where P^ is expressed in 
MPa*» the resultind r is expressed in cm/s* 
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5.1.4 Performance Parameters 

The theoretical performance parameters of the motor based 
upon the two particle flow modelsr as well as the actual results 
are presented in Table 5.2. These results are based upon the 
motor eouipped with the standard conical nozzle and the 65/35 0/F 
ratio propellant. 



parameter 


THEORETICAL 
PARTICLE MODEL 


THEORETICAL 
PARTICLE MODEL 

Tjtsconsl. 


ACTUAL 

VALUE 


THRUST (MAX), N, 


1697 


1153 


1155 








0 378 


Vj , mys. 


1490 


1821 






1406 


1406 


1347 


P^,MPa 


10.55 


10.SS 


^ " 


Cf 


2.36 


2.00 


— 


1 , N-s. 


40S 


226 


28B 


>sp, 


1 


103 


131* 


X 


0.9 89 


0.969 






0.998 


0.998 




tv 






0.633 



















Table 5.2. Comparison of theoretical performance 

parameters to actual performance parameters. 

The enerSy conversion efficiency is based upon the actual 
performance of the motor as compared to the maximum performance 
possibles which would be attained for the first particle flow 



model 
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CHAPTER 6J DISCUSSION OF RESULTS 
6*1 Nozzle Testinsi 

Contrary to expectations » the performance of the cubic 
nozzle was somewhat lesser to that of the conical nozzler deliv- 
ering a specific impulse 2*5 percent lower* This result however» 
is tentative for only a sinsle firing of the cubic nozzle was 
performed due to time limitations* Clearly^ however? no signifi- 
cant increase in nozzle efficiency is implied* This sud^ests 
that the conical nozzle is certainly a satisfactory desian. Con- 
siderina the far Greater simplicity of fabrication? the conical 
desian would have to be considered the superior choice for the 
amateur rocket enthusiast* 

If further testina reinforces the present findinas? the 
lesser performance of the cubic nozzle could be a result of its 
areater maximum diveraence anale (20 dearees v*s* 12) and/or its 
shorter lenath <17% shorter)* The latter of these would be con- 
sistent with the particle flow theory where a decrease in perfor- 
mance would be a direct result of particle laa* In a shorter 
nozzle the particles would have less time to accelerate to the 
aas velocity before exitina the nozzle* This would suaaest that 
for a small motor with a hiah mass fraction of particles? the 
nozzle should have a small diveraence anale in order to increase 
the lenath. The existina diveraence anale should possibly be re- 
duced to iO or 8 dearees. As well? the effects of reducina the 
converaence anale in order to increase the overall lenath of the 
nozzle should be investiaated * Possibly the best approach to the 
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desisin of 3 nozzle for such conditions would be to shape the con- 
tour in a manner such that the maximum flow acceleration is re- 
ducedjp to minimize particle laS. 

6*2 Motor Testinsi of VeryinS 0/F Ratio Grains 

Maximum performance for the hollow cylindrical erain is 
achieved with a moderately fuel rich 0/F ratio* TestinS indicates 
that the specific impulse is highest in the range of 65 to 70 
percent oxidizerr tapering off with a lower oxidizer percentage^ 
and being drastically lowered at high oxidizer ratios* The sen- 
sitivity at high oxidizer ratios is shown by the results that a 
maximum impuse is achieved at 70 percent oxidizerr but by incr- 
easing the ratio by 5 percent results in an immeasurable impulse* 
The highly fuel rich ratios <0/F < 55/45>also resulted in 
greatly reduced performance* despite theoretical predictions that 
the specific impulse would remain essentially constant over the 
entire range from 40 to 70 percent oxidizer* This could be attri- 
buted to the observation that a large fraction of the solid matter 
(mostly carbon) remained in the combustion chamber in the form of 
a solid porous mass* This is precisely the conditions of the 

second particle flow model where v ss constant* resulting in 

s 

3 specific impulse reduced by a factor of < 1-X ). 

At very high 0/F ratios (0/F > 70/30)* the reduction in 
performance could as well be explained by this particle flow model* 
Examination of the exhaust blast area revealed that the potassium 
carbonate particles ejected from the nozzle were actually in the 
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form of sizeable droplets* This would create the conditions of 

V R} constantr where the relatively larste particle size would not 
s 

allow acceleration of the particles to a hislh velocity* The for- 
niation of the droplets could be explained thusly* durins combus- 
tion a voluminous Quantity of potassium carbonate is formed> beins 
in the liauid state* As flow occurs through the nozzle multiple 
collisions occur between the liouid particles^ resulting in co- 
hesion and effective Srowind in particle size until either the 
maximum stable size is achievedr or exit occurs from the nozzle* 
The practical ranslef therefore? for 0/F ratio selection 
appears to be an oxidizer percentage of between 55 and 70 percent* 
Coupling this result with the realization that actual combustion 
temperatures appear to parallel theoretical? would indicate that 
since lower combustion temperatures misht be more desirable? the 
lower 0/F ratios misht be more suitable* If however? the combus- 
tion temperatures encountered at the upper end of the 0/F ratio 
pose no practical problem? the rande for best performance could 
be reduced to the 65 to 70 percent oxidizer rande* Since no 
erosion or other problems have been encountered due to these tem- 
peratures? this would appear to be the optimal rande* One final 
factor to be considered is the castability of the drain* Since 
a more sucrose rich propellant is easier to cast? the best over- 
all ratio would then appear to be the existind 65/35 0/F ratio* 

It should be noted that since testind occurred in steps of 5 
percentade points? that the optimum ratio is therefore within 
a certain plus-minus rande of this value* 
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6»3 Burnrate Testing 

The results of testing indicate that burnrate is Quite 
sensitive to the 0/F ratioi* with the maximum occurring at the 
45/35 0/F ratio* For small rocket motorsf a greater rate of bur- 
ning is desirable to develop the high chamber pressures that such 
motors can readilu operate at* With lower burnratesr it would be 
necessary to increase the grain surface (burning) area to compen- 
sate for the reduced rate of burning in order to maintain the gas 
generation rate* A decrease in throat area would as well be a 
solution to developing the desired pressure* however* a corres- 
ponding thrust loss would result* 

Increased grain surface area can be achieved bu modifying 
the geometry of the central bore* fay moulding a star or other 
shape as discussed in section 2*2*3* This result can also be 
achieved by increasing the length/diameter ratio of the grain* 
However* both methods complicate the grain casting procedure and 
increase the likelihood of grain casting flaws* Therefore* practi- 
cal considerations make the higher burnrate that can be achieved 
with the 65/35 0/F ratio a desirable trait* 

The method of estimating the burnrate-pressure relation- 
ship appears to be satisfactory for a first approximation. However* 
it is difficult to gauge the accuracy at this point* Knowledge of 
the actual chamber pressure would certainly increase the precision 
of this technioue* Due to the nonlinearity of the eauations* the 
value for the chamber pressure* however* could be subject to app- 



reciable error 
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6*A Performance Parameters 

From table 5.2 f it can be seen that there is an appreci- 
able difference in ideal thrust for the two particle flow models. 

In the second case where the particles are assumed to have no app- 
reciable velocity» the thrust is 32 percent lower than the first 
case where all the particles are assumed to have a velocity eoual 
to that of the Sas flow. As a result of this rather significant 
variancer the choice of model will importantly affect the expected 
performance of a siven motor desi:Sn» Neither model would be ex- 
pected to accurately represent the actual flow conditions since 
these are the two limiting cases. The actual particle flow cond- 
itions would reauire a consolidaton of the two modelsr where the 
expected performance would lie somewhere between these two bounds. 
As wellr the actual flow would be difficult to model since the 
real particle size would certainly be distributed, where the smal- 
ler particles mi^ht follow closely the eas flow velocity whereas 
the larger particles mi^ht have a significant velocity Iasi. 

The actual thrust produced by the motor would therefore 
be expected to lie somewhere between these two ideal values, which 
it in fact does. The actual maximum thrust value <1155 N.) almost 
coincides with the ideal value for the second particle flow model 
(1153 N.)f suSSestinS that the actual flow would be accurately 
described by this model. However, such a conclusion would be hisihly 
tentative, is not incorrect, since the calculated thrust values 
for these two models are almost directly proportional to chamber 



pressure. The chamber pressure used for these two calculations 
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(lOfSS Mpa«) is based on the value arrived at through the pressure- 
burnrate calculations described in the previous section whichf as 
mentioned* could be subject to appreciable error. An« such error 
would be almost directly reflected in the ideal thrust compu- 
tations. This factor renders the thrust parameter a poor basis of 
comparison of the actual to theoretical flow conditions. A far 
better comparison would be made utilizing the specific impulse* 
which is insensitive to chamber pressure variations. 

The value of the actual average specific impulse (131 s.) 
is in fact close to the midpoint of the values predicted by the 
two particle flow models (152. first models 103t second model). 

This implies that a consolidation of the two models would closely 
predict the behaviour of the actual flow. 

The energy efficiency e* and the velocity coefficient T * 

V 

are employed to relate the actual thrust to the ideal thrust that 
would be achieved* as explained in section 3.2.7 . The values ob- 
tained appear to be on the low side* adain possibly as a result 
of error in the chamber pressure assumed. 

The actual value for the thrust coefficient C * would 

f 

also likely lie between the two values as shown. This result 
demonstrates that the thrust can be more than doubled with a 
properly desidned nozzle* with the importance of the diverdent 
section apparent. 

The nozzle exit velocity v *is shown to be sidnif icant ly 

e 

hidher for the second particle flow model* where this represents 
the das only value. The exit velocity for the first model where 
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the particles are assumed to have the same velocity as the das 
flow is reduced as a result of acceleratind the particles* 

The actual combustion temperature <1347 C.) was found 
to be in close adreement with the theoretical value <1406 C«) beinsf 
4 percent lower. The slidhtly lower combustion temperature would 
be a result of incomplete combustion and nonadiabatic conditions. 

The final parameters to be considered are the discharde 

coefficient f. and the diverdence correction X. It can be seen 

d 

that only minor thrust losses are a result of these two factors. 
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CHAPTER 7X CONCLUSIONS 

Theoretical analysis of a small solid pro»»ellant rocket 
motor was conductedn with the performance results presented in 
the form of two limiting models* The actual findings were ex- 
pected to lie within this randef and experimental measurements 
found this to be the case* This result demonstrated the impor- 
tance of considering the effects of two-phase flow for small 
solid propellant rocket motors. 

Theoretical analysis suddested that the propellant per- 
formance should remain essentially constant over a fairly wide 
ranae of oxidizer/fuel ratios? however? other factors such as burn- 
rate and castability reduced the acceptable rande* The best over- 
all ratio was found to be the presently employed ratio of 65/35* 
Sidnificant improvement of the motor's performance does 
not appear likely to be achieved by variation of the nozzle de- 
sign* The conical nozzle emerded as a highly satisfactory design 
when the simplicity of fabrication is considered* Marginal in- 
crease in nozzle efficieny might be obtained by reducing the con- 
vergence and/or divergence angles* 

A thrust-time curve of the rocket motor's performance was 
obtained? providing a basis upon which to determine the expected 



acceleration? velocity? and altitude of the rocket in actual flight* 
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RECOMMENDATIONS 

The followinsl recommendations are suddested to provide a 
more complete understanding of the rocket motor's and propellant's 
performance? and to provide a sound basis for the design of 
larger rocket motors* 

1* Conduct testing in order to determine the actual cham- 
ber pressure durinsf firinsir possibly utilizing a pressure trans- 
ducer system to obtain a pressure-time curve. As well? a method 
of determining burnrate as a function of chamber pressure should 
be investigated utilizing this data. 

2. Other methods of determining the pressure-burnrate rel- 
ationship should be considered? such as construction of a constant 
pressure vessel to study the rate of burning at elevated pressures 

3. Investigation into modifying the conical nozzle design 
in order to reduce particle acceleration might be attempted. 

4. Further studying of two-phase flow in an attempt to con 
solidate the particle flow models would certainly prove useful. 

As well? it is suggested that a computer program be de- 
veloped to solve the eauations of motion involved with a racket 
vehicle in flight? utilizing the obtained thrust-time data in con- 
junction with vehicle drag data (either experimentally obtained or 



theoretically modeled) 
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NOMENCLATURE 

a burnrate coefficient 

a acoustic velocity 

A area 

A burning area 
b 

A/D analogue / digital 
AFT adiabatic flame temperature 
♦ 

A critical throat area 
c effective exhaust velocity 
C. Celsius 

C Coefficient of thrust 
f 

C Specific heat at constant pressure 

p 

C Specific heat of solid (lieuid) 
s 

C Specific heat at constant volume 

V 

CV characteristic velocity 

D srain inside <bore) diameter 
i 

D arain outside diameter 
o 

e eneraw conversion efficiency 
f mass fraction 
F thrust 
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a acceleration of gravity 
h enthalpy 

h“ enthalpy of formation 
f 

I total impulse 
I specific impulse 

SP 

k ratio of specific heats 

K Kelvin 

L grain length 

m mass flourate 

m gas mass flowrate 

g 

m particle (solid) mass flowrate 
s 

m.Pt melting point 

M Mach number 

M^ molecular weight 

n burnrate exponent 

N Newton 

0/F oxidizer/fuel 

P pressure 

P ambient pressure 

a 

P nozzle exit pressure 
e 

P stagnation pressure 
o 

r burnrate 

R molar gas constant 




universal sias constant 



time 

burntime 

temperature 

stagnation temperature 
velocity 

nozzle exit velocity 
sas velocity 

particle (solid) velocity 
chamber volume 

particle (solid) mass fraction 
diversience andle 
discharge coefficient 

velocity coefficient 

diversfence correction factor 

density 

Sas density 

propellant density 
stagnation density 



micrometre 
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APPENDIX A 

Calculation of the Adiabatic Flame Temperature assuming 
the combustion eouation shown belowt 



+ 3»14 K^COg +3.14 Nj 



with the assumptions of no work? no changes in K.E. or P.E.f 
and adiabatic conditionsr the first law reduces to. 



H = H 
R P 

R P ® 

where the enthalpies of formation are Siven below. 



Constituent 


State 


h| (kJ/kmol) 


C ^2 ^ 22 % 

KNOg 


solid 

solid 


“2,222 ,100 
~ 493,205 


CO 


das 


“ 110,529 


COg 


das 


- 393,522 


a 

CM 

X 


das 


“ 241,827 


KgCOg 


liouid 


“1 , 146,835 


^2 


^3S 


0 



Reference 

C313 

C323 

C33J 

C333 

C333 

C32D 





The values for Ah are zero for the reactants. For the 
productsr the values f or A h are available from standard enthalpy 
tables (except for K^COg ) such as from reference C333. The value 
for KgCOg is obtained from the expression. 
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where C = 

p 

and Ah = 
tr 



Ah 



■[> 



C dT + h 
p t r 



0.0948 T + 94.25 kJ/kmol-K C323. 



27,633 kJ/kmol C323, 



This leads to» for K2C0g. 

Ah = 1^(0.0474 T + 94.25 T) ~ 32,29 a] + 27,633 kJ/kmol 

Insertion of above values and the expression for enthalpy 
for the potassium carbonate sives the following expressionJ 

2 

8.29 < Ah ) + 0.564 (Ah > + 11 ( Ah ) + 3.14 (0.0474 T 
CO COg HgO 

+ 94.25 T) + 3.14 (Ah ) - 2,099,714 = 0 

A solution for T (AFT) is obtained throuSh a trial and 
error method by inserting values for enthalpies at the trial 
temperature. 

Solution of the eouation is obtained at T ~ 1681 K. 
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APPENIiIX B 

Determination of the averasie value for the ratio of 

specific heats kr for the 65/35 0/F ratio for flow through the nozzle. 

The assumption is made where* P = 10*55 MPa* 

o 

P - 0*101 MPa* 
e 



T = 1629 K. 
o 



Use of the expression* 



n. 


c + 


! 

1 


c + 


i 

C i 
1 


c + * * * 












% 



and the expression for C as a function of temperature siven 

p 

below C3T. 

0*75 0*5 -0*75 

COJ C = 69*145 - 0*7046 0 ~ 200*77 0 + 176*76 0 

P 

0*5 2 

CO^: C - -3*736 + 30,529 0 - 4*1034 0 + 0*024198 0 



0*25 0*5 

H-Q: C = 143*05 - 183*54 0 + 82,751 0 

P 



3.6989 0 



-1*5 -2 -3 

N, * C = 39*060 - 512*79 0 + 1072*7 0 - 820*4 0 

P 



where 0 = ? valid in the rande of 300 - 3500 K 

100 



k can be expressed as k = 
das constant* 



f where R' is the universal 



Cp- R‘ 



The Sbove expression for C becomes* 

p 



p 23 



<8*29 C + 0*564 C + 11 C + 3*14 C ) 



P, 



CO 



CO, 



HjO 
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The use of the isentropic relation! 




allows the temperature of the Sas flow at nozzle exit to be detei — 

mined» The specific heat is calculated at T and T and an average 

, o e 

value determined! 



avsi. C = ( C + C ) /2 
P P T P T 

o e 

whereas the average value for k can be found from the expression 

relating C to k> Howeveri* since it 'is the average value of k 
P 

that is to be used in the isentropic expression to obtain a true 

value for T t succesive iteration is reouired to obtain the de- 
e 

sired average value for k for flow through the nozzle! 



■^0 


Te 


avsat. Cp 


BBDI 


1420 




41.97 


1.247 


1420 


444 


38.05 


1.280 


1420 


588 


37.74 


1.282 


1420 


583 


37.74 


1.283 


1420 


583 







The value of T has reached a constant valuer so the 
e 

value for average k for flow throuah the nozzle is! 



k =• 1.283 
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APPENBIX C 



Derivation of the expression for v for two-phase flow 

e 

for the case where T « T and v rs v ♦ 

s s S3 

Consider a frozen flow model where particles are formed 

within the combustion chamber with no subseouent drowth or form- 

♦ ♦ 

ation* Under such conditionsf m and m are constant throusihout 

s M 

the nozzle* Applying the sleneral momentum relation to an incre- 
mental control volume Sives the result* 



-AdP =m <v +dv)+m(v +dv>-mv -mv <1) 
ss s ss ss dsi 

where A is the cross sectional area at a Siven point. From 

continuity. « . 

m = p A V and m = p A v (2) 

s s s S S S 

Combining these two eouations we obtain* 

-dP sspvdv + pvdv <3) 

s s s S d si 

The enerSy eouation for steady isentropic flow relates 
the particle and saseous enthalpies and kinetic energies accor- 
ding to* 

♦ ♦ 

m <C dT + V dv ) + m <C dT + v dv > = 0 (4) 

sss ss sP3 as 



Usinsi the mass fraction of particles X» this eouation can 
be re-written as* 



X C dT 



+ (1-X) C dT + X V dv 
P d 



+ <1-X) V dv = 



<5) 



s s 



s s 



0 
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From eausiion <3) we obtaini 

dp P, dP X 

V dv “ - ~ V dv = “ V dv (6) 

a sf Fg s s J-X a a 



Substitution of this expression into <5) aives* 



X X dP 

C dT + C dT + (v ” V ) dv = (7) 

1-Xss Pa 1-Xsa s Fg 



The solution of this eoustion reouires knowledae of the 

heat transfer and draa processes in order to express T and v 

5 s 

in terms of the overall flow properties* Howeveri* if the heat 

transfer is considered to be very fast (T « T )f and particle 

s a 

draa very hiah <v « v )» then a limitina solution can be obtained* 

s a 

The eauation would then reduce to* 

X dP 

< C + C ) dT = --- (8) 

1-X s P "g 

assumina the ideal aas relation? o ~ P/RT? this becomes* 

'a 



[X <l“X)] Cg + Cp 

R 



dT 

T 



dP 

P 



<9> 



Intearation from staanation conditions to exit conditions 



aives* 




where m 



R 



[X <1”X>] C + C 
s P 



For these, conditions? eauation <5) becomes 
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fx C + <1-X> C 1 dT + V dv = 0 

■- S pJ 



which may be intearated to aive v in terms of T and T t and where 

e o e 

C = k R / <k-l) : 

P 




1 “ 




m 



where 



m 



X C 

<1-X) R' 



+ 



k 

k-1 






^1 
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APPENDIX D 

Derivation for the expression for thrust F» for two Phase 

flow where T « T and v w v 4 
s d s a 

The expression for das mass flow rate throusfh a nozzle is 
aiven by the expression C34I{ 




For the case where all particles are accelerated throudh 
the nozzle (none remainina in combustion chamber) f the expression 
for total mass flowrate m» is aiven byJ 

* * 4 

m = m + Bi 
a s 



mj X 

since -4 — = » this would aive 

m 1-X 

a 



♦ 

m 



4 4 X 

m + m 

a a l“X 



4 1 

m 

a i“X 



Eauation (1) can be used to express the total mass flowrate 
in terms of the aas wherei 
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Since F 

expression for v 




where m = 



* m V + <P - P )Ar usinS the previously derived 

e 

e 6 3 

» (appendix Of we can write* 
e 




X C M' k 

+ 

(l-X)R' k-1 



-1 




E-1 



APPENDIX E 



Rocket Motor Design and Specifications 




1 3 0 

L = 18«6 cffl« L » 0.65 cm. D = 3.495 cm. 

2 4 i 



Chamber wall thicknessf t = 1.5 mm 

Chamber wall material - mild steel tubindi* s == 400 MPa. 

u 

Safetw Bolt - 3/16 in. x 1 1/2 in.r arade 6 

tensile strendthf s = 900 MPa. (rated) 

u 

shear strengths s = 448 MPa. (tested) 

s 

Nominal chamber pressure reauired to sever safety boltf R' 



p' s 2 



Nominal 



s 




where D. is bolt 

b 

chamber pressure real red to burst caseir 



diameter . 
P" 



p / y 




35 MPa. 



Dj where S is chamber tensile strength. 

u 

Nominal chamber operating pressure ~ 10.55 MPa. 
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Conical Nozzle Oesistn 




Dimensions in cm> t 



A 


3.50 


E 


0*08 


1 


0*35 


B 


0.932 


F 


2*2 


a 


12’ 


C 


2.68 


G 


1*2 


0 


30* 


D 


6*8 


H 


0*70 







Expansion Ratio! 8*27 X 1 

Material: mild steel 



Inside flow surface* polished 
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APPENDIX F 

STRAIN GAUGE BRIDGE ~ AMPLIFIER CIRCUIT 







